Waiting time distributions for the transport through a quantum dot tunnel coupled to 

one normal and one superconducting lead 
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We have studied the Waiting Time Distributions (WTDs) for sub-gap transport through a single- 
level quantum dot tunnel coupled to one normal and one superconducting lead. The WTDs reveal 
the internal dynamics of the system in particular the coherent transfer of Cooper pairs between the 
dot and the superconductor. The WTDs exhibit oscillations that can be directly associated to the 
coherent oscillation between the empty and doubly occupied dot. The oscillation frequency is equal 
to the energy splitting between the Andreev bound states. These effects are more pronounced when 
the empty state and double occupied state are in resonance. 
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Introduction — Electron transport in mesoscopic con- 
ductors is an inherently stochastic process [1] and in or- 
der to fully characterize it, it is necessary to study the 
statistics of the transport events [2] . A well established 
theoretical tool is the Full Counting Statistics (FCS) IS- 
IS]. The FCS is usually defined in the long-time limit, 
when the time interval during which transport events are 
counted is long enough that many particles have passed 
through the system. The long-time FCS describes fully 
zero-frequency transport quantities, such as the average 
current, the zero-frequency noise and higher-order cur- 
rent cumulants. Only recently the FCS has been ex- 
tended to the finite frequency domain [10-13]. Another 
tool that has recently been employed to characterize 
mesoscopic transport is the distribution of delay times 
between subsequent transport events, the Waiting Time 
Distribution (WTD) ]14-19]. When memory effects can 
be neglected, the WTDs can be used as a theoretical 
tool to evaluate the zero-frequency FCS ]14, 17]. How- 
ever, their utility lay in the fact that they are particularly 
suited to study the short-time behaviour of the system. 
This is particularly important for systems that have in- 
ternal dynamics. For this reason WTDs have been em- 
ployed to study double-quantum dots [14-16, 19], where 
coherent oscillations between states localized in the dif- 
ferent dots can occur. Quantum dots contacted with su- 
perconducting leads [20, 21] can show an interesting dy- 
namics due to the coherent exchange of Cooper pairs be- 
tween the dot and the superconductors. In this particular 
case the coherent oscillations occur between dot states 
with different particle numbers. The microscopic mecha- 
nism underlying this effect is Andreev reflection ]22, 23], 
which leads to the appearance of sub-gap resonances in 
the density of states, so called Andreev bound states, 
which have been measured by means of transport spec- 
troscopy ]24; 25]. In a quantum dot with Coulomb repul- 
sion that is tunnel-coupled strongly to a superconducting 
lead and weakly to a normal lead, a finite pair-amplitude 
can be induced in the dot, which is facilitated by the 
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Figure 1. (a) Schematic setup: quantum dot (QD) tunnel- 
coupled to one normal lead (N) with tunnel-coupling strength 
Fjv and one super conducting lead with tunnel-coupling 
strength Fg. Super conducting lead and QD are considered 
as a combined hybrid system. The normal lead is model as an 
electron bath, (b) Time-line with jumps processes, a jump Jo 
to the single-occupation sector always has to be followed by 
a jump Jo to the even-occupation sector. The general wait- 
ing time w{t) is the distribution of all Aim, m £ N. In the 
exemplary sequence of jumps shown, the waiting time matrix 
entries We,o(t) and Wo,e(i) are the distribution of all At2m-i 
and At2m- The two curves in time intervals Aii and At4, are 
sketches of the WTDs. 



non-equilibrium due to finite applied voltages [26-28]. 
Such a pair amplitude, describes the coherent exchange 
of Cooper pairs between the dot and the superconductor 
and should also be visible in the WTDs. In the present 
letter, we calculate the WTDs for a single-level quantum 
dot tunnel-coupled to one normal and one superconduct- 
ing lead in the uni-directional transport regime, see Fig- 
ure 1(a). We find that the WTDs show oscillations which 
are a signature of the coherent transfer of Cooper pairs 
back and forth between dot and superconductor. Finally, 
we discuss possible ways to measure the WTDs without 
destroying the coherent oscillations in the quantum dot. 
Model — Since we are interested in the sub-gap An- 



dreev transport between the QD and the superconductor 
rather than in quasi-particle transport, we assume the 
superconducting gap to be large. In this hmit, the dot 
in proximity to the superconductor is described by the 
effective Hamiltonian [28, 29] 



Hcs^s^dldc 



+ Un^ni - —{d 



\dl + d^d^), (1) 



with d\{da) the creation (annihilation) operator of an 
electron with spin a =t, ^ in the dot and ricr the corre- 
sponding number operator. The single particle energy 
e is assumed to be independent of the spin and U de- 
notes the on-site Coulomb repulsion. The last term in the 
Hamiltonian Eq. (1) is due to the coupling to the super- 
conducting lead and accounts for the coherent tunneling 
of Cooper pairs in and out of the dot. The prefactor Vs is 
the coupling strength between the dot and the supercon- 
ducting lead. The eigenstates of i?eff with an odd number 
of electrons are the singly-occupied dot states \a) = fi|^|0) 
with energies e, where |0) is the empty-dot state. Those 
with an even numbe r of elec trons are t he Andre ev-bound 
states |±) = -1=0^^10) T ^^^l±^\D), with 

eigenenergies e± = 6/2±eA^ where \D) ~ dldAfS). Here, 
5 ^ 2£ -\- U \s the detuning between empty and dou- 
ble occupied state, and 2eA = y/^"^ + Tg is the splitting 
between |-|-) and |— ). The normal lead is described by 

^N = Y.ka^NkC^ka'^Nkc,, where c\i^^[cNka) is the cre- 
ation (annihilation) operator of an electron with quan- 
tum number k and energy e^k- The electron distribu- 
tion in the normal lead is the Fermi function /(w) = 
(^Qp(<^-p.) _|_ \)-^ where /i is the chemical potential and 
/3 = {ksT)'^, with T being the temperature and kg the 
Boltzman's constant. The proximised dot is coupled to 
the normal lead by means of the tunneling Hamiltonian 
Ht — '^ka ^wcjvfecr'^cr + H.C., whcrc Vn is the tunneling 
amplitude. We define the tunneling coupling strength 
Tn = 27r:^N|^Np, where vj^ is the density of states in the 
normal lead, which is assumed to be constant. We trace 
out the reservoir degrees of freedom of the normal lead to 
obtain a generalized master equation [8, 27, 30] for the re- 
duced density matrix of the system, p. Since the waiting 
times are defined for unidirectional transport, we set the 
chemical potential of the normal lead to be much larger 
than all relevant energy scale of the system (apart from 
the superconducting gap). In first order in the coupling 
strength Fn the generalised master equation is Marko- 
vian and in Liouville space it can be written [6, 10] as 
p = Wp, where W is the Liouvillian or kernel. 

Introducing a counting field x [3, 4, 6, 8, 12] which 
counts the electrons tunneling out of the normal lead 
into the proximised dot, we can write the counting- 
field dependent equations of motion as p{x) = (Wo + 
e~^^J')p{x), where J is the jump operator describing 
the particles tunneling into the system from the normal 
lead. The kernel Wq describes the dynamics of the sys- 




Figure 2. The waiting time distributions We,o{T), «;o,e(T), and 
■w{t) for 5 = OU, OAU, lU. For ease of visualization, Wb,o{t), 
Wo,c{t) have been multiplied by a factor 1/2. For comparison, 
the dashed curves for Wo,e{T) and it)(r) are computed without 
the inclusion of the off-diagonal elements of the reduced den- 
sity matrix between the Andreev bound states. Parameters: 
Tn = O.Olf/ and Ts = 0.2U. 



tern in interaction with the bath without the jumps and 
contains both terms in zeroth and first order in Fn. The 
explicit expression for the kernel Wo and the jump oper- 
ator ^7 can be found in the Supplementary information. 
Results — The waiting time distribution, that is the 
probability distribution, for two subsequent jumps to be 
separated by the time interval r when the system is ini- 
tially in the stationary state described by ps, is given by 
the expression[14-16, 19]: 



w{t) 



TTJJe'^o^Jps 
nJps] 



(2) 



The system under consideration is a multi-reset sys- 
tem, since it can host more than a single excess electron. 
For multi-reset system, in order to be able to recover the 
long-time FCS one identifies different types of jumps and 
defines a matrix W{t) [14] with elements Wk,i- 



Wk,l{T) 



Tr[Jfce^""J;Ps] 
Tr[J,p5] 



(3) 



and J = ^j, Jk- The entries of the waiting time matrix 
are WTDs for a jump of the type Jk given that the last 
jump was of the type Ji. For the hybrid system under 
consideration we define two different jumps: Jo which de- 
scribes jumps from the Andreev bound states |±) (even- 
occupation sector) to the single-occupied states \a) (odd- 
occupation sector) and Jc describing the opposites pro- 
cess. An example of sequence of jumps for our system is 



shown in panel b) of Fig. 1: notice that the jumps Jo and 
J7e alternate regularly. The waiting time matrix in our 
case has the dimension 2. Instead, of working in the time 
domain we find more convenient to express our results 
in Laplace space; the Laplace transform of an arbitrary 
function gij) is defined as g{z) = J^ dre^^'^glr). Since 
in our system, subsequent jumps of the same type are im- 
possible, we find Wo, 0(2) = Wo,oiz) = 0. The off-diagonal 
elements are given by 



Wo,o(z) 
We,o(z) 



TnTI (Tn + z) 



T%rl + z^a2 + 2zrNai + Az^Tn + z* ' 



r 



N 



(4a) 



(4b) 



with ai= S^ + rjj+ r^, a2 = 6^ + SF^ + F^. Here, it 
is important to emphasize that in order to describe cor- 
rectly the short-time dynamics of the system, it is neces- 
sary to include the off-diagonal elements of the reduced 
density matrix between the Andreev bound states, which 
we refer as the coherences. Since in the hybrid system a 
jump to the single-particle sector |ct) has to be followed by 
a jump to the |±) sector and vice versa, the general wait- 
ing time of Eq. (2) reads w{t) = i (we^oiT) + Wo,c{t))- 

The waiting time Wc^o{t) is given by the expression 
Wc,o{t) = Fjv exp (— Fatt). This can be easily understood 
since after a jump jTo the dot is in a singly-occupied state 
and therefore decoupled from the superconductor. The 
next jump will then be a Poissonian process with rate F^r. 
The situation is different for Wo,e(T). The jump J will 
bring the system from \a) to \D) which is a coherent su- 
perposition of the two eigenstates |-|-) and |— ). Therefore 
we expect coherent oscillations with a frequency that in 
zeroth-order in Tn is given by the Andreev-bound state 
splitting 2eA- Although we have the full result in Eq. (4), 
it is still instructive to look at the two limiting cases when 
the proximity effect is in resonance, 6 = 0, and off reso- 
nance, (5 ::^ F5. On resonance, 2eA = Fg and the waiting 
time Ws^aIt) reduces to 



.(r) 



_ SnT^-f^t 



3=0 



F2 



F2 



1 



cos 



r| 



N' 



Fjve-^"ni-cos(F5r) 



(5) 



where in the second equality we have made use of the 
fact that Ftv ^ F5. On the other hand, off-resonance, 
2£a ~ 1^1; we find for short times tTn <C 1 



= (t) 



TatF 



5, 



-tFj, 



(52 



[1 - cos(|<5|r)] 



and for long times tFjv S> 1 



W^o,o(t) 



JVi^^-. 



2(52 



(6a) 



(6b) 



The off-resonance long time behaviour is consistent with 
the picture of Poissonian tunneling of Cooper pairs found 




Figure 3. Time evolution of Pdd + poo, Pdd and poo obtained 
by the time evolution with the Liouvillian without jumps Wo 
when the system is initially prepared in the doubly-occupied 
state for 5 = QU , O.lf/, If/. Parameters: Fat = O.Olf/ and 
Fs = 0.217. 



in Ref. [28]. The full result for the waiting times w{t), 
Wc.oiT) and Wo,ciT) obtained by the inverse Laplace 
transform of Eq. (4) is shown in Figure 2 for different 
detunings 5, including intermediate ones. The WTD 
'Wo,c{t) always start from zero for t = 0, since after the 
jump J7c the system is in \D) and any further transport 
is blocked until the two electrons in the dot are trans- 
ferred coherently to the superconductor. We notice that, 
as the detuning increases, the amplitude of the oscilla- 
tions is suppressed and the frequency of the oscillations 
is increased. From the discussion above it is clear that in 
order to describe the oscillations it is necessary to include 
the coherences between the Andreev-bound states. Not 
including the coherences gives only the exponential be- 
haviour and does not resolve the oscillating terms. This 
is demonstrated in Figure 2, where we also show w(t) and 
Wo.c{t) obtained without the inclusion of the coherences. 

In order to understand the features of Wo,ciT) we in- 



vestigated the time evolution of p{t) 



oWot 



PD , that 



is the evolution of the system without jumps when at 
t = the system is in the doubly-occupied state. The 
quantity p{t) are the elements of the density matrix 
under the condition that no jumps have occurred in 
the time interval t. We use the localised basis for the 
dot and write p = {poo, Pn^ PU^ Pdd, Pdo, Pod}, with 
PTt(^) ~ Plli't) = since no jump has occurred. The 
time evolution of pdd and poo clearly governs the evo- 



lution of Wo^dT). Figure 3 shows the time evolution of 
Pdd + Poo, Pdd and poo for different detunings S. As the 
system is initially in the doubly occupied state we have as 
initial values pDD{t=0) + poo(i=0) = 1, PDD{t=0) = 1 
and poo{t=Q) = 0. All the elements of p are decaying 
since the probability of having no jumps decays with time 
and eventually goes to zero as times goes to infinity [31]. 
On resonance the local maxima of pooit) are equal to 
PDD{t) + poo{t) at all times, since the coherent dynam- 
ics govern the evolution of p and the behaviour of the 
WTDs at all times. This is not the case off resonance in 
the long-time limit {tT^ S> 1), where the coherences are 
suppressed and the WTDs show only a Poissonian decay, 
see Eq. (6b). 

The waiting time matrix in Laplace space W{z) can be 
used to derive the long-time FCS. In particular, the cu- 
mulant generating function is the solution of the equation 
[14] 



0.020 



det[e^x - W{z)] = 0, 



(7) 



t hat fulfils zo(x=0) = 0. We find zq = -Fat + 

j= , which re- 
duces in the limit of small Fat to the result of Ref. [28] 
obtained without the inclusion of the coherences between 



|±), zo 



N 



S^i 



-2^xTl 



s^+n 



- 1 



We also define the conditional waiting time as 



wcA^)^TT[Je^«-p^] 



(8) 



which is the probability distribution for the next jump 
to happen after the time r given that the system has 
been prepared initially in the state |^). The condi- 
tional waiting time can be used to measure the entries 
of the waiting-time matrix directly. Since after a jump 
from the odd sector to the |±) sector is in the double- 
occupied state, we have that the conditional waiting time 
u!c,\D){z) is simply wc,\D)iz) = Wo^dz)- Similarly, we 
find wc.\a-) (z) = Wc.oiz), with the difference that a single 
occupied state has no coherent evolution in time, such 
that only the next jump changes the system state. The 
conditional waiting time distribution for the system ini- 
tially in the empty state is an independent quantity from 
the entries of the waiting-time matrix and reads 



Wc,\o)iz) = 



Tn {{Tn + z) {2z (F^ + z) + F|) + 25^z) 



^ N'^ S 



z^a2 + 2zFArai + Az^T 



N 



Figure 4 shows the conditional waiting times for the sys- 
tem being prepared in the empty state Wc^\o)iT) and the 
double-occupied state wc,\D)i''')- Obviously u'c,|d)(0) = 
since the system is in the double-occupied state and 
transport is blocked, while 'Wc,\o)i^) is at the maximum 
value since the system is in the empty state and by con- 
struction no electron could have tunneled at an earlier 




Figure 4. Conditional waiting times Wc,|o> i^) (system initially 
empty |0)) and wc,\D)iz) (system initially double-occupied 
\D)) for (5=0(7 and S^OAU at riv=20rs. The dashed curves 
are calculated without the inclusion of the coherent dynamics. 



time as T = 0. We also show the curves obtained with- 
out inclusion of coherent dynamics. 

Finally, we discuss the possibility to measure the 
WTDs for the system under consideration. The WTDs 
for uni-directional transport though a quantum dot can 
can be measured by monitoring the occupation of the dot 
by means of a charge detector, e.g. a nearby capacitively- 
coupled quantum-point-contact [32-37[. In the present 
system, the situation is slightly more complex as Cooper 
pairs oscillate coherently back and forth between the 
dot and the superconducting lead and a charge detector 
would decohere these oscillations. The ideal detector to 
measure the WTDs in this system is a device that mea- 
sures whether an unpaired spin is present on the dot. As 
a model let us assume that the current in the detector is 
of the form Idet ~ isin-^ + n^^'2n^n^). The detector cur- 
rent is equal to ig if the dot is in the single occupied state 
and it vanishes if the dot is in the even-occupation sector. 
The detector current as a function of time corresponds to 
a time trace for the jumps Jc and j7o. The distribution 
of the length of the time intervals with current is is the 
WTD Wc.o (r) : similarly the distribution of the length of 
the time intervals with zero current is the WTD Wo,c (''')• 
Since such a detector is not able to discriminate between 
the states |0) and \D), it does not introduce decoherence 
in the coherent tunneling of Cooper pairs. 

Conclusion — In this Letter we have calculated the 
WTDs of a single-level quantum dot coupled to one 
normal and one superconducting lead in the regime of 
strong-coupling to the superconductor. The WTDs re- 
veal features which are directly related to the Andreev- 
bound states. In particularly, we found that the WTDs 
oscillate with a frequency equal to the Andreev bound 
state splitting. The amplitude of the oscillation is maxi- 
mal when the proximity effect is on resonance and probes 
the pair amplitude in the dot. The time-scales associated 



to the coherent Cooper pair oscillations do not appear in 
the zero-frequency FCS. The entries of the waiting-time 
matrix can be measured directly with a device that senses 
the presence of an unpaired spin in the dot. 
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SUPPLEMENTARY INFORMATION 

The generalised master equation for a nanostructure tunnel-coupled to a fermonic bath is in general non Markovian 
[1, 2]. It becomes Markovian when the coupling to the bath Tn = 27r^j. \Vn\'^S{uj — SNk) is taken into account in 
first order in perturbation theory or in the high-bias limit (/i — ;■ ±00) when electronic transport is unidirectional. 
Therefore, for the regime considered here, the kernel is Markovian and the generalised master equation is of Linblad 
form and can be written as 



p = Wp=-i[Hs,p]+V{p), 



(S.l) 



where Hs is the system Hamiltonian, ^{p) the Lindblad dissipator and we have set ft = 1. For the system under 
consideration, we have Hs = -ffoff and in the high bias limit {p — > 00) the dissipator reads 



T^ip) 



Y, ^N (e~'^dlpd^ - hpd^di + dadlp) 



(S.2) 



where we have introduced the counting field x that counts electrons tunneling out of the normal lead. The first part 
of the dissipator describes a measured change in the number of particles of the normal lead, i.e. a jump event. The 
second part describes dissipation of the system due to the coupling to the normal lead without jumps. The kernel is 



then rewritten as W(x) = Wq 



^J, where J' is the jump operator. We separate the jump operator J = Je + Jo, 



into a part Je describing jumps to the the Andreev-bound states (even) sector and a part Jo describing jumps to the 
single-occupied (odd) sector. To make this separation we introduce a different counting field for each type of jump 
and we rewrite the kernel as 



W{xe,Xo) = Wo + e-^>^= Je + e-'^^-X. 



(S.3) 



The kernel takes a simple form in the basis {poo, Ptti Pll^ Pdd, Pdo, Pod}, where p^/^ ~ (C'IpIO ^^d the states |^) and 
^') belong to the localised dot basis {|0), | t), I i), \D)}. With this choice of basis in Liouville space, the kernel reads 
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(S.4) 



From Eq (S.4), Wo and the jump operators Je and Jo can be read off. 

For the sake of completeness we also give the expression for the kernel in the basis of the eigenstates of Hcg: 
{I t), I i), l+)l^)}- Using the basis {p^i^i, Pn, p ,/?++,/? y,P-i } in Liouville space, we obtain 
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This representation of the kernel can also be easily obtained by means of a diagrammatic real-time expansion [3] . 
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